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ABSTRACT 

Satellites  are  implicitly  completely  described  by  their  orbital  elements.  The 
conditions  for  going  from  the  mean  orbital  elements  of  a  satellite  to  osculating 
elements  considering  perturbations  caused  by  the  ellipsoidal  earth,  and  then  to 
celestial  coordinates  and  their  rates  of  change  are  derived.  For  belts,  it  is 
necessary  to  fix  a  point  at  which  it  is  desired  to  direct  an  antenna.  This  is  done 
by  taking  the  intersection  of  a  right -ascension  half -plane,  a  geodetic  longitude 
half-plane,  a  declination  cone,  or  a  time  varying  central  angle  vector  with  the  belt. 
The  equations  for  determining  this  point  from  the  mean  orbital  elements  and  the 
intersecting  surface  or  ray  are  derived. 
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PREFACE 

This  document  is  a  collection  of  notes  by  W.  G.  Williams  (Parke  Mathematical 
Laboratories,  Carlisle,  Mass.)  that  were  used  by  C.  W.  Adams  Associates,  575 
Technology  Square,  Cambridge,  Massachusetts  in  constructing  two  programs  for 
the  Haystack  Pointing  System  under  a  subcontract  with  Group  62  of  Lincoln  Laboratory. 
The  document  was  prepared  by  C.  W.  Adams  Associates  and  is  a  quite  literal  trans¬ 
cription  of  handwritten  notes.  Included  also  is  an  appendix  of  editorial  notes 
amplifying  certain  sections  of  the  notes  or  noting  de  facto  changes.  These  sections 
are  indicated  by  an  asterisk. 
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INITIAL  TRACKING  PROGRAM  CALCULATIONS 
(June  25,  1963) 


Preliminary  Remarks 


A.  The  coordinate  system  which  Smithsonian  Astronomical 
Observatory  (SAO)  employs  is  the  customary  sidereal 
system  referred  to  the  equinox  of  1950  and  the  equator 
as  of  the  date.  The  coordinate  system  employed  by 
Lincoln  Laboratory  is  the  customary  sidereal  system 
but  referred  to  the  equinox  of  date  and  the  equator  of 
date.  In  the  SAO  system  the  measurement  of  the  right 
ascension  of  the  ascending  node  is  from  the  equinox  of 
1950  along  the  fixed  equator  to  the  ascending  node  of 
the  moving  equator  and  thence  along  this  equator  to  the 
ascending  node  of  the  satellite  (Kozai,  A.J.,  ^5,  621, 
1960)  .  No  problems  should  be  encountered  by  making 
only  a  simple  correction  to  the  ascending  node  of  the 
satellite  to  take  into  account  the  precession  of  the 
equinox  alone  (that  is,  negating  nutational  effects), 
since,  to  within  high  precision,  both  equators  (SAO 
and  LL)  will  coincide  after  the  precession  term  has 
been  built  in. 


B.  The  SAO  mean  elements  will  be  given  as  'a  function  of 
time,  generally  having  the  form 


and  occasionally  a  term  of  the  form  D  sin  m  or  E  cos 
to  take  into  account  the  effect  of  the  third  harmonic  of 
the  terrestial  gravitational  field. 

Time  t  in  this  expression  is  given  by 


t  =  T  -  T0 


and  is  in  days,  where  To  is  the  reference  epoch  (that  is, 
a  certain  Julian  Day  and  fraction  thereof)  in  days. 
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C.  The  outline  of  the  problem  is  as  follows: 

1.  Once  having  been  given  the  date  and  the  particu¬ 
lar  satellite  to  be  followed,  evaluate  the  mean 
elements  (whose  form  has  been  cited  in  B)  for  an 
initial  time  T. 

2.  From  this  evaluation  calculate  the  short  period 
perturbations  (which  will  be  given)  and  add  to 
these  mean  elements  to  get  the  instantaneous  os¬ 
culating  elements  at  the  time  T. 

3.  From  these  osculating  elements  calculate  the  right 
ascension  and  declination  and  radius  of  satellite. 

4.  Proceed  to  the  time  T+&T  and  repeat  this  process. 

Note:  It  may  turn  out  that  it  will  be  sufficient  to  cal¬ 
culate  the  osculating  elements  only  once  during 
the  given  pass.  If  so,  then  item  2  above  will  be 
calculated  this  once  and  then  we  will  assume  that 
bo  in  the  general  expression  of  B  is  this  oscula¬ 
ting  element  and  that  b  =  b  (t)  will  be  the  time 
change  to  a  sufficient  precision  for  the  pass. 


D.  The  elements  employed  by  SAO  are  the  usual  elliptic  set: 
uu  =  argument  of  perigee 

n  =  longitude  (right  ascension)  of  ascending  node 
i  =  inclination 
e  =  eccentricity 
M  =  mean  anomaly 

n  =  anomalistic  (mean)  motion  =  M 
a  =  semi-major  axis 

Units  are :  degrees  for  angular  quantities 

mega-meters  (106  meters)  for  linear  quantities 
revolutions  for  M  and  derivatives  of  M 


2 


N.B, 


SAO  employs  a  Modified  Julian  Date  which  is  related  to 
the  Julian  Date  as  follows:  MJD  =  JD  -  2,400,000.5 


Calculation  of  Satellite  Right  Ascension  fop  . 

Declination  (6)  and  Radius  fp). 

Given  the  SAO  mean  elements  (as  described  above)  and  the 
initial  time  of  the  calculation  (i.e.,  T  which  is  a  J.D.  and 
some  fraction)  and  the  reference  epoch  (which  may  be  an  M.J.D. 
and  some  fraction)  to : 

1 .  Calculate  MJD  of  T  and  then  calculate  t  =  T  -  To  . 

2 .  Evaluate  the  orbital  elements : 

ou  =  U)  (t) 

n  =  n(t) 
i  =  i(t) 
e  =  e(t) 

M  =  M(t) 
n  =  M(t) 

3.  Correct  0  for  the  precession  of  the  equinox  by  add¬ 
ing  to  Q  the  term  3.S08xl0~5  (MJD-33281),  i.e., 

“  5 

n  =  n  +  3.508x10  (MJD-33281) 

4.  Evaluate: 


(GM  is  a  constant  which  is  to  be  specified) 
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5  „  Evaluate : 


a 


which  can  be  approximated  as : 


a 


3  .  s  . . 

2  sm  t) 


(As  is  a  constant  to  be  specified) 

6.  Calculate  the  eccentric  anomaly  E  from  the  Kepler 
equation : 

E  =  M  +  e  sin  E 

This  will  be  done  by  an  iterative  process ;  as  an 
initial  trial  value  use  the  expression: 

3 

E  =  M  +  (e  -  g-  )sin  M  +  ^e2  sin  2M  +  §e3  sin  3M  +  0(e4) 

7 .  Calculate  the  true  anomaly  v  from  the  relations 


cos  E-e 

1-e  cos  E 


sin  E;  cos  v  = 


8 .  Calculate  argument  of  latitude  u  from  the  expression 
u  =  uo  +  v 


9  o  Evaluate  p  =  a  (1  -e  cos  E) 


10.  Calculate  short  period  perturbations  from  the  ex¬ 
pressions 

dfi  =  -  cos  i  jV-M+  e  sin  v  -  isin  2  (v-ku) - 

^  sin  (v+2uu)  -  ^  sin  (3v  +  2uu)j- 

1  As  f 

d£  =  ■jj  2  sin  21  jcos  2(v-K«)+e  cos  (v+2tu)  + 

P 

•j  cos  (3v+2uu)j- 

du  =  p|*  [^(2-  |  sin2i)  (v-M+  e  sin  v)  +  (1-  sin2t) 

(1-  -  a/i -e2  sin  v  +  |r  (1-  Vl-e2  sin  2vj- 

-  C|-  -  ^  sin2t)  e  sin  (v+2uu)  ~  “  f2  sin2t) 

sin  2  (v+w)  -  ^  cos2  i  sin  (3v+2uu)  J 
dp  =  (1-  sin2t)  j^-1  -  (1-  ,/l-e2)  cos  v+ 

+  f  -Tu?  }  +  ?  ^  slnH  cos  2<v4<"> 

11 .  Apply  perturbation  in  range  (from  center  of  earth) 
p  =  a(l-e  cos  E)  +  dp 

12.  Apply  perturbation  to  node, 

fi  -*  +  dH 
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13 „  Apply  perturbation  to  inclination  in  the  form 

sin  i  sin  z  +  cos  i  dz  1  where  trig  functions  on 

.  .  .  .  |  right  are  unperturbed, 

cos  z  -*  cos  z  -  sm  z  dz  J 

1M-.  Apply  perturbation  to  argument  of  latitude  in  the  form 

sin  u  sin  u  +  cos  u  dzz  ]  ,  .  .  -  .  . 

I  where  trig  functions  on 

.  ,  right  are  unperturbed 

cos  u  -»  cos  u  -  sin  u  du  ■ 

J 

15.  Calculate  the  right  ascension  (a)  and  the  declination 
(6)  from  the  relations 

sin  6  =  sin  z  sin  u  =  Ab 

cos  a  cos  8  =  cos  fi  cos  u  -  cos  z  sin  0  sin  u  =  Ax 

sin  a  cos  ft  +  sin  ft  cos  u  +  cos  z  cos  ft  sin  u  =  Ab 

Calculation  of  Pass  Visibility  at  This  Time  T 

1  .  Calculate  the  right  ascension  of  Greenwich  at  the  time 
T  from  the  Newcomb  Formula 

round  off  to  10  places 
i  l 

aG  =  [0*2769193981+0*002737909263  (T-2,415,020)  + 

+0  .0000010752  (3r~6525xl0^0)3]  ^  ^in  radians) 

*  2  .  Calculate  right  ascension  of  radar  a  at  time  T  from 
aR  ‘  aG  +  AR 

where  =  longitude  of  radar  (West  is  negative ; 

East  is  positive) 
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*  3 .  Calculate  declination  of  radar  from  the  relation 

2 

6  =  t  -  m  sin  2t  +  - —  sin  4£ 

K  d 

where  -t  =  geographic  latitude  (North  taken  positive) 
m  =  3  .359076852xl0“3 

4 .  Calculate 


Bi 

=  COS 

6r 

COS 

a 

R 

Bb 

=  cos 

6r 

sin 

a 

R 

% 

=  sin 

6r 

5.  Calculate  T,  the  angle  between  satellite  direction 

(from  earth  center)  and  radar  direction  (from  earth 

center)  from 
3 

cos  r  =  E  A.  B. 

•  a  7/7/ 


6 .  Convert  geographic  minimum  "look"  angle  3  at  radar 

to  geocentric  "look"  angle  3  by  means  of 
/  2 
3  =  3  +  m  sin  2t  -  —  sin  4-t 


*  7  .  Calculate  geocentric  distance  of  radar  from  the  re¬ 
lation 


R  =  Rp  [1  +  (K  cos2£)  +|  (3C  cos2t)2]  +  h 

where : 


R 

P 

h 

K 


=  polar  radius  of  earth 

=  height  of  radar  above  sea  level 
=  3  ,34783xl0”3 
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*  8  .  Calculate  visibility  zone  half  angle  G  from  the 
relation 


cos  G  =  f~)  cos2 8  +  sin(3^  \l  1  -  (jr)2  cos2 3^ 

where  r  is  the  geocentric  radius  of  satellite . 

*  9  .  Satellite  is  visible  when 
T  *  G 
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INITIAL  BELT  TRACKING  PROGRAM 
(September  9,  1963) 


The  following  theoretical  treatment  is  based  upon  the 
fact  that  the  same  tracking  accuracies  will  ultimately  be 
required  of  the  belt  program  as  are  required  of  the  satel¬ 
lite  program;  and,  further,  that  the  orbit  parameters  will 
in  effect  be  given  in  terms  of  SAO  mean  elements  .  If  the 
second  assumption  is  unfounded,  the  analysis  will  be  simpli¬ 
fied  considerably  in  an  obvious  way  since  a  large  share  of 
the  analysis  consists  of  going  from  mean  to  ’’instantaneous" 
osculating  elements  in  a  straightforward  way.  If  the  first 
assumption  can  be  realized,  the  number  of  iterations  re¬ 
quired  is  obviously  less . 


*  Determination  of  First  and  Last  Belt  Visibility  Times 

Oth  Approximation 

1 .  Calculate  cos  G  (given  in  satellite  Initial  Track¬ 
ing  Program,  which  shall  be  referred  to  as  ITP) 
from  equation  8  of  the  Visibility  section  with 
r=a  where  a  is  gotten  at  5  of  the  first  section 

of  ITP.  Convert  to  sin  G,  noting  that 
0  £  G  ^  tt/2. 

2.  Calculate  the  mean  elements  of  the  belt  at  the  be¬ 
ginning  of  the  J.D.  (Julian  Day)  of  interest. 

3  .  Determine  the  first  and  last  visibility  right  as¬ 
censions  of  the  radar,  a  ,  from  the  equalities  of 
the  following  relation (in  which  6^  is  the  declina¬ 
tion  of  the  radar  calculated  in  3  of  the  visibility 
section  of  ITP) : 
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-sin  G  ^[sin  i  cos  6  sin  (D-a  )  +  cos  i  sin  6D] 

K  K  K 

^  sin  G . 

Note  that  there  will  be,  in  general,  two  such  sets 

of  contacts  per  sidereal  day  because  of  the  rotation 

of  the  earth;  thus  when  the  two  a^'s  have  been  found, 

there  will  be  two  more  derived  by  adding  fi  +  •?  -  a 

<-  K. 

to  fl  +  Tf  in  the  Oth  approximation . 

4 .  The  universal  times  for  these  contacts  can  be  de¬ 
rived  from  these  a  ' s  by  solving  for  a  in  ITP  visi- 

K  '  b 

bility  2,  and  then  inverting  the  Newcomb  Formula  of 
ITP  visibility  1 . 

This  completes  the  Oth  approximation 

5.  To  start  the  1st  approximation,  solve  the  follow¬ 
ing  expression  for  the  argument  of  the  latitude,  u, 

of  the  belt  at  these  contact  a  ’  s  from  the  expres- 

K 

sion : 

cos  6  cos  (a  -Q)  cos  u  +  [cos  6  cos  i  sin  (a.  -Q) 

K  K  K  K 

+  sin  6  sin  i]  sin  u  =  cos  G 


1st  Approximation  Series 

6.  From  the  T’s  in  4  update  all  the  mean  elements  to 
the  appropriate  times  . 
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7 .  At  each  of  these  times  use  the  mean  elements  updated 
in  6  and  the  u's  calculated  in  5  and  calculate  the 
geocentric  distance  of  the  belt,  p,  from  the  formula: 

p  =  —.afl-e3).- 

M  1+e  cos  (u- ou) 

8.  With  these  values  of  p  recalculate  cos  G  (and, 
hence,  sin  G)  by  replacing  a  in  the  Oth  order  expres¬ 
sion  by  p  . 

9 .  Repeat  steps  3 ,  4  and  5  deriving  new  contact  right 
ascensions ,  times  and  arguments  of  latitudes . 

10 .  Repeat  the  cycle  from  6  through  9  until  consistency 
is  attained. 

2nd  Approximation 

11  .  From  the  mean  elements  finally  settled  upon  in  10 

and  from  the  relation  that  v  =  u  -  uu  where  u  is  got¬ 
ten  from  10  and  from  the  relations  that  p  = 
a  (1-e  cos  E)  and  M  =  E  -  e  sin  E,  calculate  the  values 
of  v  and  M  (which  are  needed  in  6Q  and  6 u)  and  calcu¬ 
late  the  short  period  perturbations  of  ITP  (namely, 
dfl,  di ,  du  and  dp)  .  Apply  (i.e.,  add  as  in  ITP) 
these  mean  elements  settled  on  in  10  thereby  deter¬ 
mining  the  osculating  elements  . 

12.  Repeat  steps  3,  4  and  5,  thereby  determining  the 

contact  right  ascensions  of  radar,  the  times  of  con¬ 
tact  (these  are  the  crucial  times)  and  the  argu¬ 
ments  of  latitude  of  belt  at  contact . 
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13 .  Using  these  osculating  elements  of  11 ,  determine 
the  contact  right  ascension-declination  sets  of 
the  belt  at  contact  from  15  of  ITP. 


Determination  of  Right  Ascension-Declination  Values  of  Belt 

inside  Visibility  Times 

For  a  given  visibility  interval  the  procedure  is  as 

follows : 

14.  Jump  time  from  1st  contact  by  At  and  update  mean 
elements  to  this  time . 

15  .  Solve  5  for  "upper"  and  "lower"  argument  of  lati¬ 
tude  of  belt  using  the  value  of  p  in  cos  G  derived 
from  11 . 

16.  Apply  short -period  perturbations  to  these  "upper" 
and  "lower"  points  to  get  osculating  elements  . 

17.  Using  these  values  of  the  elements,  again  solve  5 
for  the  "end  point  u's". 

18.  Using  these  u's  solve  for  the  end  point  right 
ascensions  and  declinations  by  means  of  15  from  ITP. 

19.  Using  "appropriately"  averaged  elements  (i .e . ,  the 
mean  of  the  end  point  osculating  elements)  calculate 
the  right  ascension  and  declination  of  the  belt 
(and  radar  range  if  desired)  for  a  range  of  u's 
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intermediate  to  these  end  point  u's  from  15  of 
ITP  (and  from  7  above)  . 

20 .  Jump  time  from  1st  contact  by  2At  and  cycle  from 
15  through  19. 

21.  Jump  to  3At,  4At . . . .etc . ,  until  end  of  visibility 
interval . 


Remarks 


a)  The  right  ascension-declination  values  derived  in 
19,  20  and  21  are  sufficient  to  construct  a  pre¬ 
cise  schedule  for  the 

1)  fixed  latitude  (declination)  schedule 

2)  belt  scan  schedule 

3)  fixed  right  ascension  schedule 

4)  fixed  longitude  schedule 

For  example,  a  precise  plot  of  the  values  in  a  form 
such  as  a  vs .  6  with  T  as  a  parameter  contains  all 
the  information  necessary. 

b)  The  above  treatment  has  been  done  in  complete 
generality.  For  certain  "degenerate"  situations 
and/or  with  reduced  accuracy  requirements,  it  can 
be  considerably  simplified  with  little  effort  if 
desired . 
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ADDITIONS  TO  INITIAL  BELT  TRACKING 
PROGRAM  AND  MISCELLANEOUS  ADDENDA 
(September  10,  1963) 


Initial  Belt  Tracking  Program 


*  1 .  If  instantaneous  osculating  elements  are  employed  as 

inputs  rather  than  SAO  mean  elements,  the  calculations 
are  considerably  simplified  in  that  the  process  is 
ended  at  10  at  the  termination  of  what  has  been  desig¬ 
nated  the  1st  approximation  series.  Further,  in  steps 
1-10  replace  "mean  elements"  with  "osculating  elements" 
mentioned,  and  use  "osculating"  for  "mean"  in  14  and 
delete  16  and  17  which  are  rendered  unnecessary. 


2.  In  terms  of  the  desired  schedules  the  calculation  of  18 
is  superfluous.  All  that  is  required  is  the  calculation 
of  the  "end  point  uTs".  Thus  the  output  of  the  program 
will  be  as  follows : 


T 

u 

Initial  contact 
time 

Ti 

Ui 

Intermediate 

Ti  =Ti  +AT 

u  l 

Ui  ,  Uj. 

times 

=Ti  +2AT 

u  1 

^3  5  Ua 

• 

• 

Final  contact 
times 

T 

f 

Uf 

• 

Initial  argument 
of  latitude 


Upper  and  lower 
(  u +  1)  u' s  at  in¬ 
termediate  times  . 


Final  argument  of 
latitude . 
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The  interval  AT  should  be  chosen  such  that  a  smooth 
and  precise  interpolation  scheme  can  be  worked  out  for 
times  and  u's  in  between  these  calculated.  Thus,  schemat¬ 
ically  we  may  represent  the  output  in  the  following  graph¬ 
ical  form: 


From  these  visibility  "footprints"  or,  more  precisely,  the 
appropriately  interpolated  u  vs.  T  outline  of  the  "prints", 
the  belt  schedules  can  be  worked  out  as  follows : 


A.  Fixed  Latitude  (Geodetic-)  or  Declination  Schedule 

Given:  6 to  remain  fixed 

u 

a)  from  sin6  =  sini  sinu  ^ - (I) 

solve  for  Ug  corresponding  to 

b)  From  "foot  print"  solve  for  visibility  times,  i.e.. 


t 


visibility  true  for  6  constant 
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c)  From  either 

eoscc  cos6  =  cosfi  cos u  -  cost  sinfi  sinu  ^ 

or  sinoc  cos6  =  sinfi  cosu  -  cost  cosfi  sinu 

with  6=6  and  u  =  u 
15  Jd 

solve  for  a  . 

Jd 

d)  From 

p  =  a(l-ef) _ 

1+ecos  (u- ou)  "<1 — 

with  u  =  U-  solve  for  pn 
Jd  Jd 

*  B .  Constant  Belt  Scan  Schedule 

Given :  ^  =  k  to  remain  fixed  where  v  =  true  anomaly 

rj 7/  /}  •  • 

(Note:  (v-hju)  =  k  +  uu  where  uu  =  perigeal  pre- 

•  • 

cession;  thus  u  can  be  calculated  with  k  +  w  given  as 

they  are.  However,  what  is  probably  meant  by  belt  scan 

d n  ° 

is  that  —  =  k.  At  any  rate,  one  can  get  to  u  from  k 

in  a  straight  forward  fashion .) 

a)  At  T  =  Ti  we  have  from  which  (at  ,  6t  ,  pt)  can  be 
gotten  by  means  of  (I)  ,  (II)  and/or  (III)  and  (IV)  . 

b)  At  T  =.Ti  +  AT  we  have  U  =  +(^)  At  from  which  the 

next  (a,  6,  p)  set  can  be  calculated. 


(II) 

(III) 


(IV) 
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*  c)  Repeat  Ob)  until  belt  is  no  longer  visible  at  horizon 

*  d)  Reverse  sign  of  and,  repeating  (b)  scan  "backwards” 

to  other  horizon 

e)  Continue  c)  and  d)  until  time  has  run  out. 

(This  process  is  pictured  below) 


T.  T„  T  - 

1  f 

C .  Fixed  Right  Ascension  Schedule 

Given  a  •-  to  remain  fixed 

JD 

(Note:  When  0^  =  0,  ,  TT,  3n/2  ,  u  assumes  certain  de¬ 

generate  values ;  this  situation  must  be  handled  separate¬ 
ly  ,  obviously) 

a)  Set  a  -*  and  from  the  following  (III  divided  by  II) 

tan  a  =  cosit  -  cost  cosO  sinu 

an  cosQ  cosu  -  cost-  sinQ  sinu  ^ _ (V) 

calculate  u  -*  uD 

D 

b)  Solve  for  visibility  times  as  in  A. b) above 

c)  Solve  for  6^  and  from  (I)  and  (IV) 


17 


D .  Fixed  Longitude  Schedule 


Given  Ot^  =  +  A  «-  Fixed  longitude  (W  =  -  ;  E  =  +) 

T 

Greenwich 

which  varies  as  a  function  of  time,  being 
given  by  the  Newcomb  Formula 

At  each  time  step  repeat  a) ,  b)  and  c)  of  C,  the 
fixed  right  ascension  schedule. 


General  Notes 

I)  Let  X  be  angle  between  orbit  plane  (osculating)  and 
meridian 

then 

cos  X  =  sint  cos  (ct^  -  fi) 

where  =  right  ascension  of  meridian  of  interest  (such 
as  the  ctg's  of  above)  . 

II)  It  should  be  emphasized  that  the  Belt  Tracking  Program 
provides  a  beautifully  simple  way  of  calculating  satellite 
visibility  times  as  follows: 

a)  Calculate  the  visibility  "footprint"  for  the  orbit 
of  the  satellite  thought  of  as  a  belt,  in  precisely 
the  same  way  (including  short  period  perturbations) 
as  previously  given. 
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b)  Now,  noting  that  u  =  mod  2tt,  calculate  u  as  a 
function  of  time  from  the  start  of  the  J  .D .  of  in¬ 
terest  to  end  by  means  of  the  formulas  of  satel¬ 
lite  tracking  program . 

c)  Compare  these  with  "footprints",  when  u  falls  in¬ 
side  satellite  is  visible,  when  outside  invisible. 
This  can  be  represented  schematically  on  next  page . 
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Notes : 


1.  Radar  assumed  in  Northern  hemisphere 

2.  Orbit  assumed  unperturbed 

3.  B  is  mirror  image  of  A 

w.r.t.  J.D.  +  \  true 

4.  Perigee  assumed  at  m  =  5 

5.  8  hr.  orbit  assumed 

6.  A  and  B  not  drawn  to  scale 


SATELLITE  VISIBILITY 
PERIODS 


T 


PRE-DOPPLER  CALCULATIONS 
(February  6,  1964) 

Satellite/Belt  Doppler  Effect 


Assumptions 

We  will  neglect  ionospheric,  tropospheric  and  relativ¬ 
istic  considerations : 

a)  For  further  discussion  of  these  see: 

Mass  and  Vassy;  Chapter  1  of  Advances  in  Space 
Science  and  Technology,  4,  1-38  (1962) . 

b)  For  a  more  detailed  analysis  of  the  first  and 
second  items  see: 

1)  Grier  and  Weiffenback,  Proc .  I.R.E.  48 , 

507-516  (1960) 

2)  von  Handel  and  Hohndorf,  I.R.E.  Trans  on  Mil. 
Elec.  3_,  162-172  (1959) 

c)  Note  that  neglect  of  relativistic  effects  is 

V  2 

~(— )  in  relative  frequency  shift  and  hence  will 
not  become  "appreciable”  (~  1  cps)  until  fre¬ 
quencies  of  l-lOGC(KMc)  are  employed  with  the  ve¬ 
locities  (~  iOkm/sec)  "normally"  encountered  in 
satellite  technology . 
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Theory 


It  can  be  shown  that  the  1st  order  Doppler  shift  is 
given  by  the  following  expression: 


(1) 


in  which 

Af  .=.  Doppler  shift  in  frequency 
f  .=.  Frequency  employed  by  transmitter 

Note :  It  is  assumed  here  that  there  is  an  active  transmitter 

on  board  the  satellite .  If  we  are  dealing  with  a  passive 
satellite  or  with  belts,  i.e.,  bouncing  e-m  radiation  off 
the  satellite/belt ,  then. 


o 


c 


-(Af/pc  =  Pi  +  p2  where  px 

P2 


=.  transmitter/belt  range 
=.  belt/receiver  range 


velocity  of  light  (vacuum) 


c 


range  (from  transmitter  to  receiver  if  trans¬ 
mitter  on  satellite  or  (see  note  above)  range 
from  radar  to  belt  or  from  belt  to  receiver) 


P 


o 


range  rate  in  which  a  dot  denotes  time  differ¬ 
entiation  . 


P 


Now ,  if 


r  .=.  vector  pointing  from  center  of  earth  to  satellite 
(or  to  some  specific  position  on  belt) 

— ♦ 

R  .=.  vector  pointing  from  center  of  earth  to  radar 
transmitter  (or  receiver) 
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r  .=.  angle  between  these  two  vectors  (Note:  this 
is  the  same  angle  r  specified  earlier) 

then  we  see  that  with 

cosP  =  sin6sin6R  +  cosScos6Rcos  (a-aR)  (2) 

where  the  6's  and  a’s  have  the  usual  meaning  of  declinations 
and  right  ascensions,  respectively;  that 

p2  =  r3  +  R2  -  2rRcosr  (3) 

where  cosr  is  given  by  (2) . 

Performing,  now,  the  indicated  time  derivatives  we  see  that 
P  =  ■“  |^[r-Rcosr]  r  -  rR  |^[sinSR  cosS-cos5R  sin&  cos(a-aR)]6 

[cosbR  cos&  sin(a-aR)]  (a-aR)J  j  (4) 

in  which  we  have  noted  that: 

i  =  SR  ,  0 

i  e  ,  the  "radar"  remains  on  a  small  circle  R  away  from  the 

center  of  the  earth  at  declination  6  moving  only  due  to  the 

earth's  rotation  aD . 

K 

•  •  •  • 

We  need,  therefore,  r,  &,  a  and  aR  in  addition  to  tha  param¬ 
eters,  all  of  which  have  been  calculated  before  except 
for  p  which  is  calculated  from  (3) 


We  note  that 

•  a(l-e2)esinv  * 
r  =  — - - - - -  v 

(1+ecosv)2 


(5a) 
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where 


’Vl-e2  • 


V  = 


r2 


M 


and 


c*  sinZ/  cosu  .• 

6  =  - t -  u 

coso 


where 


u=—  (U)  +  v)  =  UJ  +  V 


and 


*  A  _  A  .  cosZ/ 

CL  —  { l  +  - p-r-  H 

COS  o 


(5b) 


(6a) 


(6b) 


(7) 


which  can  be  gotten  from  Teoste ' s  expression  after  very 
tedious  algebra  or  by  noting  that : 


a  =  Q  +  cos 


cosL 

cos6 


(8) 


and  finally, 

aR  =  1  .002737909263  X2tt  radians/day  (9) 

(i  .e  . ,  7  .29211585x10" 5radians/sec) 

NOTE:  Both  the  Newcomb  Formula  and,  hence,  the  following 

relation  in  the  paper  of  June  25,  1963  are  in  error. 
See  the  accompanying  correction . 


****************************** 


Hence,  it  is  seen  that  p  can  be  calculated  after  some  effort 
by  using  that  which  we  knew  before,  for  the  most  part. 

CALCULATION  PROCEDURE  FOR  DOPPLER  SHIFT  AT  TIME  T 

I  Extract  everything  from  satellite  right  ascension  and 
declination  and  radius  portion  of  Initial  Tracking  Pro¬ 
gram  (June  25,  1963);  extract  all  but  6,  8  and  9  of 
Pass  Visibility,  at  the  time  T  (using  corrected  1  and 
2,  of  course). 

II*  Decide  (in  case  of  belt)  what  mode  to  be  in: 

a)  Fixed  declination: 

JB  =  UB  s  0 

b)  Constant  belt  scan: 

•  • 

U  =  R+uu  or  =  K 

c)  Fixed  R.A.  and/or  Fixed  Longitude:  both  determine 
5^  and  Ug  as  a  function  of  time. 

*  III  Calculate  r,  6,  a  and  from  (5),  (6),  (7)  and  (9) 

*  IV  Calculate  p  from  (3) 

*  V  Calculate  p  from  (4) 

*  VI  Calculate  Af  or  Af/f  from  (1)  (for  active  satellite) 

*  VII  a)  Double  results  of  VI  if  transmitter  and  receiver 

are  at  same  site  (Passive  satellites  or  belts) 

b)  Do  an  additional  calculation  for  other  link. 


*  Recall  note  on  p . 22 . 
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APPENDIX 


Editor's  Notes 


p.  4  Section  6 

The  iterative  process 

E  =  M  +  e  sin  E 


has  been  replaced  by  the  highly  convergent  procedure 


M.  =  E.  —  e  sin  E. 
11  l 


AE. 

l 


M  -M. 

l 


1  —  e  cos  E 


i 


E 


i  +  1 


E.  +  A  E. 


l 


l 


The  initial  value  Eq  is  found  either  by  the  series  expansion  of  E  as  shown  or  by 
a  linear  prediction  based  upon  the  last  two  calculated  values  of  E. 

See  "Celestial  Mechanics"  by  Forest  Ray  Moulton,  the  MacMillan  Co.  ,  New 
York  (1959)  p.  164. 

p.  6  Section  2 

The  right  ascension  of  the  site  at  time  T  is  found  by  the  coordinate  conver¬ 
sion  program  using  stored  tables  of  apparent  sidereal  times.  This  value  is 
available  to  other  programs  of  the  pointing  system,  including  the  belt  and  satellite 
programs. 

p.  7  Section  3,  7 

The  coordinate  conversion  program  computes  during  system  initialization, 
the  declination  (geocentric  latitude)  and  the  distance  from  the  geocenter  to  the 
radar  plane. 

p.  8  Section  8,  9 

The  satellite  program  does  not  calculate  a  visibility  zone.  Whenever  called, 
it  computes  the  celestial  coordinates  of  the  satellite  for  the  time  given  the  program. 
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P-9 

The  method  of  determining  the  visible  portion  of  a  belt  given  here  is  not  used. 
Rather,  when  there  are  two  possible  points  on  a  belt  (as  may  happen  in  the  fixed 
declination  and  the  scan  schedules),  that  point  is  chosen  whose  geocentric  radius 
vector  makes  the  smaller  angle  with  the  geocentric  radius  vector  to  the  site. 

p.  14  Section  1 

SAO  mean  elements  are  not  used  by  the  belt  program.  Another  set  of  mean 
elements  determined  by  special  purpose  programs  at  Lincoln  Laboratory  is  used. 
Only  to  and  (t  in  this  set  have  first  derivatives,  i.e. 

to  =  to  +  to  t 
o 

0  =  t 

The  other  elements  are  fixed  (over  a  run), 
p.  16  Section  B 

Indeed,  belt  scan  is  ^  =  K 
p.  17  Section  B-c  and  B-d 

The  belt  program  looks  at  the  radar  elevation  of  the  point  of  the  belt  as  com¬ 
puted  by  the  coordinate  conversion  program.  Whenever  the  elevation  goes  from 
positive  to  negative,  the  sign  of  ^  is  changed. 

p.  23  Equation  5a 

This  equation  has  been  replaced  by 


r  = 


a  e  sin  v 

7TT?~ 


M 


which  can  be  gotten  from  William's  expression  after  very  tedious  algebra  or  by 
differentiating 


r  =  a(l  — e  cos  E) 
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p.  24  Equation  7 


Because  of  scaling  difficulties,  this  has  been  replaced  in  the  program  by 

'  »  •  cos  i 

a  cos  5  =  Q  cos  6  +  —  ^  u 


The  original  expression  was 

•  cos  Q  sin  u  +  cos  i  sin  r«!  cos  u  •  _  M  t  c  *  .  • 

sin  a  cos  6 

p.  25  Sections  III,  IV,  V,  VI,  VII 

The  calculation  of  range  and  range  rate  is  performed  by  the  coordinate 
conversion  program.  Doppler  is  calculated  by  the  interpolation  program. 
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